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Abstract 

We have investigated Rb adsorption on the Si(100) surface for 0.5 and 1 monolayer 
coverages using the total energy method with norm-conserving pseudopotentials. 
For 2x1 reconstruction at 1 ML coverage symmetrized dimers are found to be 
energetically more favorable. On the other hand, half a ML coverage is found to 
have symmetrical dimers only for the most stable adsorption model. All possible 
surface-adatom configurations have been considered in the calculations to find which 
adsorption sites are energetically favored. In addition to the structural properties, 
the interface is investigated electronically for the work function and surface states. 
The results are discussed and compared with the existing experimental findings. 
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1 Introduction 



Silicon surfaces covered with alkali metal (AM) atoms have been investigated 
extensively for various coverage models. The motivation behind these studies 
is to understand the metal-semiconductor junctions and to study metallization 
of semiconductor surfaces. One may expect alkali metal adsorbed silicon sur- 
faces to exhibit similar characteristics. However, systematic differences arise 
as one goes from lighter atoms to heavier ones leading to different electronic 
behavior. The discussion, especially, on the atomic and electronic structure of 
alkali metal adsorbed Si(100)-2xl surface is still under debate which concerns 
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the matters such as the adsorption sites, the saturation coverage at room tem- 
perature (RT), the nature of alkali metal-Si bond and the charge transfer from 
the adsorbate to the substrate. 



Another reason why AM adsorption on Si surfaces attracted much interest 
is that AM's have simple electronic structures and do not intermix with the 
surface Si atoms, except Li, allowing us to investigate a relatively simpler 
interface. Therefore Si-AM systems were seen as a prototype to understand 
the physics of AM overlayers on semiconductor surfaces. In the last years the 
advances in the engineering of nano scale structures allow device designs based 
upon these contacts on the atomic scale. 



There are, fundamentally, two proposed coverage models which have been 
experimentally realized. These are 0.5 and 1 ML coverages. Half a ML model 
has been first studied by Levine for Cs/Si(100)-2x 1 system [5]. In this model 
alkali metal atoms are placed in the middle of neighboring dimers in the same 
row. This adsorption location is called as pedestal site, p, (see Figure 1). Alkali 
metal atoms adsorbed on pedestal site form a linear chain along the dimer row 
preserving 2x1 surface symmetry. 



In recent years, there is much interest towards larger alkali metals (Rb and 
Cs) [1,2,3,4]. Etelaniemi et al. studied the adsorption sites of Rb adsorbates 
on Si(100) surface with XSW experiments for various coverages [6]. They at- 
tempted to define the surface geometry at the saturation coverage. Castrucci et 
al. determined AM-Si bond length and the Rb site distribution on Si(001)2x 1 
surface for low (0.19 ± 0.02) coverages using XSW [7]. Chao et al. investigated 
the surface electronic structure and work function shifts as a function of the 
Rb coverage. [8,2]. Johansson et al. made a systematic study of the electronic 
structure for Rb adsorption on Si(100)2xl surface with ARUPS and IPES by 
comparing to the results of other AM's [1]. 



In this work our aim is to study Rb adsorbed Si(100) surface with 2x1 pe- 
riodicity for 0.5 and 1 ML coverages which include experimentally realized 
adsorption sites. Main focuses of this study are to determine the atomic struc- 
ture of Si surface reconstructed with Rb coverage, to discuss the saturation 
coverage with an emphasis on the Si dangling bonds and surface work func- 
tion, and to investigate the electronic structure to identify surface states. This 
is, to the best of our knowledge, the first theoretical work that investigates 
the Rb adsorption on Si(100) surface. 
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Fig. 1. On the left : Schematics are shown for the side view of relaxed clean 
Si(100)-2xl surface. On the right : Top view of Si(100) surface with all possible 
stable adsorption sites are illustrated. (Dimers are shown as symmetrical for visual 
convenience.) 

2 Computational Details 

The calculations were performed using the plane-wave expansion approach 
within the local density approximation to density functional theory as imple- 
mented in the ABINIT code [9]. We generated Troullier-Martins pseudopo- 
tentials [10] using fhi98pp [11] package with PW92 [12] exchange-correlation 
model. We treated the true electronic valence states in the calculations. Non- 
linear core-valence corrections were taken into account for the case of Rb by 
setting the cut-off radius to 3.6 bohr. To test Rb pseudopotential, Rb crystals 
having bcc structure have been reproduced with the lattice parameter 5.49 A. 
This value is 1.6 % lower than the experimental value of 5.58 A [13]. 

We carefully tested the convergence behavior of our calculations with respect 
both to the plane-wave cut-off and to the k-point grid. A series of tests showed 
that 8x4x1 Monkhorst-Pack mesh [14] with 16 Ryd cut-off gives well con- 
verged results for our self-consistent calculations. In other words, any further 
increase in these parameters would not enhance the accuracy of the calculated 
values like lattice constants, band structures etc. The convergence criterion to 
reach the lowest energy configuration is determined by minimizing the forces 
using the Hellmann-Feynman theorem. 

For surface geometry and band structure calculations we have used a slab with 
8 layers of Si. In order to prevent an interaction between the two surfaces, 
the bulk dangling bonds are hydrogenated and a vacuum region of 10 A is 
used . In geometry optimization calculations, Si atoms which stay at the two 
bottom layers are kept frozen and all the other atoms are let relax to their 
minimum energy configurations. In order to observe if the surface symmetry 
breaks we did not impose any symmetry operations other than the identity 
which naturally incorporates. 
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We performed separate self-consistent calculations to determine surface work 
function using a symmetrical slab having a thickness equal to 10 Si layers 
and being separated from each other with a vacuum region of ~ 10 A. After 
relaxing the clean surface into minimum energy geometry we freezed 4 Si 
atoms which are located in the middle of this slab. We also repeated geometry 
optimizations using symmetrical slab for all of the coverage models as done 
in the case of hydrogenated surface. These calculations produced equivalent 
results. 



3 Results and Discussion 

Although a number of experimental studies are present for Rb adsorption on 
Si(100) surface, no first principles work has been published yet. This paper 
not only studies Rb adsorption on Si(100) theoretically, but also presents 
calculated results like Rb-Si bond length, adsorption energy, etc. which have 
not been measured experimentally. In a few papers these values are estimated 
from the theoretical studies for K/Si(001) or Cs/Si(001) systems. We tabulated 
our calculated results for structural properties in Table 1. For the surface 
energetics we used the formula as described in the reference [15] to calculate 
the adsorption energies. The adsorption energy represents the negative of the 
binding energy. 

Table 1 

Calculated parameters for the atomic structure for 0.5 ML and 1 ML coverages of 
Rb adsorption on Si(100)2xl system. The lengths, d,a mcr , d^b-dimer an d Ax, are 
all in A. (iRb-dimcr is the distance from the adsorbed Rb atom to the dimer center. 
Ax is the height of the Rb atom measeured from the dimer center. 



e 


model 


^dimer 


^Rb— dimer 


Ax 


^tilting 


E ad (eV) 


0.5 


b 


2.37 


3.05 


3.01 


10.4 


0.74 




c 


2.30 


4.12 


1.67 


5.6 


1.21 




P 


2.39 


3.04 


2.36 


8.0 


1.28 




h 


2.32 


4.42 


1.26 


2.4 


1.69 


1 


b-p 


2.44 


2.75-3.02 


2.75-2.34 




0.50 




c-h 


2.34 


4.03-4.45 


1.30-1.27 




0.68 




b-c 


2.40 


3.03-4.16 


3.04-1.67 




1.11 




c-p 


2.42 


4.20-3.01 


1.76-2.34 




1.39 




b-h 


2.48 


2.98-4.42 


2.98-1.15 




1.46 




h-p 


2.49 


4.41-3.00 


1.11-2.32 




1.72 
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Our calculated clean surface dimer length is 2.27 A. It compares well with the 
experimental value of 2.20 ± 0.05 [16]. This value becomes larger because of 
the adsorption which causes Rb to make a charge transfer to the closest dimer 
Si atom. In turn, this saturation of the dangling bonds on the surface leads 
to a stretching of the dimer. When we compare the dimer lengths between 0.5 
ML and 1 ML coverages, we see that the dimer is more stretched for 1 ML 
adsorption models. 

For 0.5 ML coverage, bridge (b), cave (c), hollow (h) and pedestal (p) sites 
are found to be stable adsorption configurations. We tested the stability of 
two more sites which are the top of one of the dimer Si atoms and the site 
on the half way from p to h, so called to shallow site. None of these two sites 
were found to be stable. Each time we put an Rb atom on these locations 
they relaxed into the nearest stable sites. Energetically the most favorable 
adsorption site is h. The total energy of this geometry is 0.41 eV lower than 
the second lowest configuration p. This surface geometry has the dimer length 
of 2.32 A with a tilting angle of 2.4°. Symmetrical dimers which correspond 
to a tilting angle of less than 5° are found only for this model. The other 
adsorption configurations - b, c and p- result in asymmetrical dimers. 




Fig. 2. On the left : Illustration for top view of h-p (or double layer) adsorption 
model on Si(100)-2xl surface. (Upper is shown bigger.) On the right : Side view of 
Rb/Si(100)-2x 1 system for 1 ML h-p configuration. 

For 1 ML coverage, we have studied all the possible adsorption configurations 
and have found bridge-cave (b-c), bridge-hollow (b-h), bridge-pedestal (b-p), 
cave-pedestal (c-p), cave-hollow (c-h) and hollow-pedestal (h-p) are the stable 
site combinations. Any combination with the top and/or shallow sites which 
are as described in the case of half ML, are energetically unstable. Based on 
an estimation for the Rb-Si bond distance Etelaniemi et.al.[6] argue that two 
nearest neighbor adsorption sites cannot be occupied at the same time. We 
found that these configurations (b-p and c-h) produce the lowest two adsorp- 
tion energies. Therefore, they are energetically the least probable adsorption 
models for 1 ML coverage. 
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The 2x1 symmetry was found to be preserved for these 1 ML adsorption cov- 
erage models as suggested by the LEED result of Chao et al. [8,2] and ARUPS 
result of Johansson et al. [1]. Among these h-p is the energetically most favor- 
able configuration which corresponds to PV model of Etelaniemi et.al.[6]. In 
our calculations h-p has a total energy of 0.52 eV less than the second lowest 
energy coverage model b-h. Rb-Si bond length for this configuration has been 
found to be 3.20 A and 3.24 A measured from Rb atoms at /i-site and p-site, 
respectively, to the nearest Si atom. In this geometry symmetrically dimer- 
ized surface is energetically favored. After the relaxation the dimer length is 
measured to be 2.49 A which is close to the experimental value of 2.3 - 2.4 
A reported by Etelaniemi et al. [6]. This adsorption model is also referred as 
the double layer model. Rb atoms adsorbed at h and p-sites are separated a 
distance 1.21 A along [100] as shown in Figure 2. This difference in the heights 
is close to one Si layer which has the value ~1.20 - 1.30 A. In this model the 
nearest-neighbor distances for the Rb atoms have been predicted by Johans- 
son et al. to be 3.84 - 4.0 A [1] which is in a very good agreement with our 
calculated values of 3.82 A along [011] and 4.0 A along [011]. 

The work function is calculated as the difference between the vacuum level 
and the Fermi energy (Ep). The vacuum level is determined from the self- 
consistent, plane-averaged potential {V av ) in the middle of the symmetrical 
slabs along the direction [100] perpendicular to the surfaces. It can be obtained 
from the Poisson's equation, 

V av (x) = -47rp av (x) (1) 



dx 2 

and is given by the relation, 

oo oo 

V av (x) = - 47r / p av (x')x'dx' + 4ttx J p av (x')dx' (2) 



where 



Pav(z) = \ip(x,y,z)\ 2 dydz. (3) 



Calculated values for the plane-averaged electrostatic potential and the work 
function corresponding to different adsorption models are shown in the Fig- 
ure 3 by setting the zero level at the Fermi energy. Figure 3(a) represents clean 
surface. The work function for Si(100)-2xl surface was reported as 4.9 eV by 
Abukawa et al.[17]. It agrees very well to our theoretical result of 4.9 eV, for 
the clean relaxed surface. Figure 3(b) corresponds to the half ML model of 
Levine [5] which is the /i-model in our case. For 1 ML coverage we presented 
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the result in Figure 3(c) which represents PV model of Etelaniemi et.al.[6] 
which is the h-p model in our calculations. The zigzag regions in the middle of 
these figures represent Si layers in the bulk. A Si atom layer sits in each valley 
and a localized charge coincides with each peak in this bulk region. Each peak 
above the Fermi level represents the position of a Rb atom. The tales of the 
electrostatic potential away from the bulk region go into the vacuum regions. 




2.20 
0.0 
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Electrostatic Potential (eV) 




(b) 



1.65 
0.0 






(C) 




Fig. 3. Planar-averaged electrostatic potential V(x) (a) for the clean relaxed geom- 
etry, (b) for h adsorption model and (c) for h-p model along [100], the direction 
perpendicular to the surface within the cell, (d) Calculated shifts in work function 
as a function of the coverage. 

The work functions were calculated as 2.20 eV and 1.65 eV for the cases of 
0.5 ML and 1 ML, respectively. When compared with the results of Johansson 
et al. [1] and Chao et al. [2], we obtained a similar behavior for the change in 
the work function as Rb coverage increases. These work function shifts as a 
function of coverage are shown in the Figure 3(d). The lowering in the work 
function is more rapid between the clean surface and 0.5 ML coverage than 
that of between 0.5 ML and 1 ML coverages. This behavior is consistent with 
the experimental observations [1,2]. Moreover, they reported that the work 
function takes a constant value after the saturation coverage. This can be 
explained by the fact that there is no further lowering in the work function by 
Rb adsorption after the saturation reached for Si (100) surface. It appears from 
the Figure 3(d) that 0.5 ML cannot be the saturation coverage for Rb/Si(100)- 
2x1 system since there is a shift in the work function for the coverages greater 
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than half a ML. We calculated the change in the work function for 1 ML to 
be 3.25 eV. The experimental result for the saturation coverage is 3.4 eV 
given by Johansson et al. [1] and Chao et al. [2]. The comparison of our 
theoretical results with the experimental findings gives a strong evidence that 
the saturation coverage occurs at about 1 ML. 

The bulk projected surface band structure is presented in Figure 4. Energy 
values are given relative to the valence band maximum (VBM) along high 
symmetry points. The direct band gap is calculated to be 0.47 eV which is 
underestimated by LDA as expected. The experimental value was given as 0.6 
eV by Johansson et al. [1]. Conduction band minimum (CBM) and valence 
band maximum (VBM) both occur at T. The overall electronic behavior re- 
veals similar characteristics to the Si(100)2x 1-AM band structure except the 
dispersion and location of surface states. We identified surface states which 
fall into the band gap. Two unoccupied and two occupied surface states are 
shown in Figure 4. 



2 - 




r j k j' r k 



Fig. 4. Bulk projected surface energy bands. The shaded region is the projected 
bulk band continuum. 

Chao et al., in their ARUPS study, reported two Rb-induced surface states CI 
and C2 for a RT saturated Si(100)-2x 1 surface with energies of -0.55 and -1.55 
eV at T, respectively [2]. Their results are in favor of the proposed double- 
layer model which leads to a semiconducting surface. The states CI and C2 
correspond to the occupied states Si and S 2 in our case with the energies -0.35 
and -0.72 eV at T. 

Johansson et al. observed a Rb-overlayer derived empty band, U 2 , with a large 
isotropic dispersion [1]. They argued that even though they observed a met- 
allization due to the large dispersion of the unoccupied band U2 at saturation 
coverage, the surface electronic behavior is expected to be semiconducting at 
1 ML with a surface band gap of about 0.6 eV as they measured. The ex- 
planation of the occupation of this empty band was given by the presence of 
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small amounts of excess charge on the surface together with the close energy 
position of U 2 at VBM. The minimum of this band was reported to be about 
0.2 eV above the VBM at T [1]. There is a good agreement with this experi- 
mental value such that we found the separation between the U 2 minimum and 
the VBM to be 0.25 eV. Although they observed a metallization at saturation 
coverage we see from the Figure 4 that Rb/Si(100)-2xl system has a semi- 
conducting electronic structure at 1 ML for the h-p-model even we see a large 
dispersion of 0.84 eV in the conduction band U 2 . 



4 Conclusion 



In the case of 0.5 and 1 ML coverages, all possible configurations were studied 
for Rb adsorption on Si(100) surface having p(2xl) symmetry. For 0.5 ML 
coverage, /i-site and for 1ML coverage, h-p model are found to be the minimum 
energy configurations. Our results suggest that Rb atoms prefer to be adsorbed 
firstly, on h-site then on p-site as the coverage increases. 

Clean surface work function was calculated to be 4.9 eV which is equivalent to 
the experimental value of Abukawa et al [17]. The work functions for 0.5 ML 
and 1 ML are 2.20 eV and 1.65 eV, respectively. A 3.25 eV shift in the work 
function was found for h-p model which agrees well with the experimentally 
observed shift of 3.4 eV for the saturation coverage at RT. Bulk-projected 
surface bands were calculated within the affinity of the gap for the case of 
h-p model. Experimentally resolved 5i and S 2 bands together with the unoc- 
cupied U 2 band are identified as surface states. 1 ML Rb adsorbed Si(100)-2xl 
surface was electronically found to be semiconducting. Being consistent with 
the experiments, the separation at T was found to be 0.25 eV between VBM 
and Rb driven band U 2 . When the change in work function and the energetics 
of the surfaces are compared, our results support the double-layer model (h-p) 
as the saturation coverage. Our theoretical results fit well to the picture drawn 
by the experimental findings. 
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